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We discuss the dynamics of linear, scalar perturbations in an almost Friedmann-Robertson- Walker 
braneworld cosmology of Randall-Sundrum type II using the 1+3 covariant approach. We derive 
a complete set of frame-independent equations for the total matter variables, and a partial set 
of equations for the non-local variables which arise from the projection of the Weyl tensor in the 
bulk. The latter equations are incomplete since there is no propagation equation for the non-local 
anisotropic stress. We supplement the equations for the total matter variables with equations for 
the independent constituents in a cold dark matter cosmology, and provide solutions in the high and 
, low-energy radiation-dominated phase under the assumption that the non-local anisotropic stress 

. vanishes. These solutions reveal the existence of new modes arising from the two additional non-local 

' degrees of freedom. Our solutions should prove useful in setting up initial conditions for numerical 

^ codes aimed at exploring the effect of braneworld corrections on the cosmic microwave background 

(CMB) power spectrum. As a first step in this direction, we derive the covariant form of the line of 
' sight solution for the CMB temperature anisotropies in braneworld cosmologies, and discuss possible 

mechanisms by which braneworld effects may remain in the low-energy universe. 



00 

(N 
> 

cn 
cn 
o 



PACS numbers: 04.50.-|-h, 98.80.Cq 



I. INTRODUCTION 



It is understood that Einstein's theory of general relativity is an efFective theory in the low-energy limit of a more 
general theory. Recent developments in theoretical physics, particularly in string theory or M-theory, have led to the 
idea that gravity is a higher-dimensional theory which would become effectively four-dimensional at lower energies, 
t Braneworlds, which were inspired by string and M-theory, provide simple, yet plausible, models of how the extra 
^ i' dimensions might affect the four-dimensional world we inhabit. There is the exciting possibility that these extra 
I , dimensions might reveal themselves through specific cosmological signatures that survive the transition to the low- 
' energy universe. It has been suggested that in the context of braneworld models the fields that govern the basic 
rr| interactions in the standard model of particle physics are confined to a 3-brane, while the gravitational field can 
^ ' propagate in 3 -t- d dimensions (the hulk) . It is not necessarily true that the extra dimensions are required to be small 
k>( ' or even compact. It was shown recently by Randall and Sundrum for the case d = 1, that gravity could be localized 
: to a single 3-brane even when the fifth dimension was infinite. As a result, the Newtonian potential is recovered on 
' large scales, but with a leading-order correction on small scales: 

-«-^(-|^)^ 

where the 5-dimensional cosmological constant A cx —l~^. As a result, general relativity is recovered in 4 dimensions 
in the static weak- field limit, with a first-order correction which is believed to be constrained by sub- millimeter 
experiments at the TeV level [|l|, ^ . 

The cosmic microwave background (CMB) currently occupies a central role in modern cosmology. It is the cleanest 
cosmological observable, providing us with a unique record of conditions along our past light cone back to the epoch 
of decoupling when the mean free path to Thomson scattering rose suddenly due to hydrogen recombination. Present 
(e.g. BOOMERANG @) and MAXIMA g) and future (e.g. MAP and PLANCK) data on the CMB anisotropies 
and large-scale structure provide extensive information on the spectrum and evolution of cosmological perturbations 
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potentially allowing us to infer the spectrum of initial perturbations in the early universe and to determine the 
standard cosmological parameters to high accuracy. An obvious question to ask is whether there are any signatures 
of extra dimensions which could be imprinted on the cosmic microwave sky. 

The aim of this paper is to set up the evolution and constraint equations for perturbations in a cold dark matter 
(CDM) branc cosmology, presenting them in such a way that they can be readily compared with the standard 
four-dimensional results, and to provide approximate solutions in the high and low-energy universe under certain 
restrictions on how the bulk reacts on the brane. Our equations are clearly incomplete since they lack a propagation 
equation for the non-local anisotropic stress that arises from projecting the bulk Weyl tensor onto the brane, and 
our solutions are only valid under the neglect of this stress. However, our presentation is such that we can easily 
include effective four-dimensional propagation equations for the non-local stress should such equations arise from a 
study of the full bulk perturbations. The lack of a four-dimensional propagation equation for the non-local stress 
means that it is currently not possible to obtain general results for the anisotropy of the CMB in braneworld models. 
Such a calculation would require solving the full five-dimensional perturbation equations which is non-trivial since 
the equations can only be reduced to two-dimensional partial differential equations on Fourier transforming the 3- 
space dependence. Only qualitative results are currently known, obtained with either the standard metric-based 
re-invariant) approach [| |, 0, |, |[ 0, |ll|, |l2|, ^ |l|, |l6|, [l7[ |l|, |o), or with 1+3 covariant methods 

In order to make this paper self-contained, we begin by giving a brief overview of the 1-1-3 covariant approach to 
cosmology and define the key variables we use to characterize the perturbations in Sec. |l[ After a short review of how 
this approach can be used to describe the general dynamics of Randall-Sundrum braneworlds, in Sec. |l^ we present 
a complete set of frame-independent, linear equations describing the evolution of the total matter variables and the 
non-local energy and momentum densities for scalar perturbations in an almost-FRW universe (with arbitrary spatial 
curvature). Many of these equations, which employ only covariantly defined, gauge- invariant variables, have simple 
Newtonian analogues p^ , and their physical meanings are considerably more transparent than those that underlie 
more standard metric-based approaches. In Sec. VII we derive analytic solutions for the scalar modes in the high 
and low-energy radiation-dominated universe neglecting the non-local anisotropic stress. In principle, these solutions 
could be used in a phenomcnological manner to generate more general solutions which include non-local anisotropic 
stress, using Green's method and an ansatz for the non-local stress. Such a study will be presented in a subsequent 



paper 1 25 , along with a numerical calculation of the CMB power spectrum employing this phenomenology. As a first 
step, we derive the covariant line of sight solution for the temperature anisotropics from scalar modes in braneworld 
models in Sec. VIII and present some comments on how higher dimensional effects may remain in the low-energy 
universe and thus imprint on the microwave sky. An appendix presents the scalar perturbation equations in the 
matter energy frame during radiation domination. 



II. 1+3 COVARIANT DECOMPOSITION 



Throughout this paper we adopt the metric signature ( h ++). Our conventions for the Riemann and Ricci 

tensor are fixed by [Va,Vb]C/c = Rabc^Ud and Rab = Racb'^- Lowercase latin indices a...b are used to denote 
the standard 4-dimensional (1+3) spacetime whereas uppercase A . . . B and the tilde of any physical quantity are 
used to denote 5-dimensional (1+4) spacetime (of the braneworld). Round (square) bracket denote symmetrization 
(antisymmetrization) on the enclosed indices. We use units with ?i = c = 1, so the 4-dimensional gravitational 
constant is related to the 4-dimensional Planck mass via G — AIp^^. 

We start by choosing a 4- velocity This must be physically defined in such a way that if the universe is exactly 
FRW, the velocity reduces to that of the fundamental observers to ensure gauge-invariance of the approach. From the 
4- velocity we construct a projection tensor hab into the space perpendicular to u°' (the instantaneous rest space 
of observers whose 4- velocity is m"): 

Kb = 9ab + UaUb , (II. 1) 

where gab is the metric of the spacetime. The operation of projecting a tensor fully with hab^ symmetrizing, and 
removing the trace on every index (to return the projected-symmetric-trace-free or PSTF part) is denoted by angle 
brackets, i.e. T(^ab...c)- 

The symmetric tensor hab is used to define a projected (spatial) covariant derviative I?° which when acting on a 
tensor T^ ''^d...e returns a tensor that is orthogonal to u° on every index, 

Z?'^r^ -^d...e = h%hK . . . h%h'd . ■ . /i"eVf T"^- , (II.2) 
where V denotes the usual covariant derivative. 
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The covariant derivative of the 4-velocity can be decomposed as follows: 

aUt = LJab + CTab + ■^'dhab - UaAb , (II. 3) 

where uJab is the vorticity which satisfies u°-ujab — 0, cr^b — 0(ab) is the shear which is PSTF, 8 = V^u^ = measures 
the volume expansion rate (where H is the local Hubble parameter), and Aa = uf'^bUa is the acceleration. 

Gauge-invariant quantities can be constructed from scalar variables by taking their projected gradients. Such 
quantities vanish in the FRW limit by construction. The comoving fractional projected gradient of the density field 
p^*' of a species i (for example, photons) is one important example of this construction: 

A« EE ^i^apW , (II.4) 

where a is a locally defined scale factor satisfying 

d EE u^Vba ^ Ha , (II.5) 

which is included to remove the effect of the expansion on the projected gradients. Another important vector variable 
is the comoving projected gradient of the expansion, 

Za = aDaQ , (II.6) 

which provides a measure of the inhomogeneity of the expansion. 

The matter stress-energy tensor Tab can be decomposed irreducibly with respect to m° as follows: 

Tab = PUaUb + 2U(^aqb) + P^ab + T^ab , (II. 7) 

where p = Tabu°'u^ is the energy density measured by an observer moving with 4-velocity u"" , qa = —h^aTbcU'^ is the 
energy flux or momentum density (orthogonal to u'^), P = habT""^ /Z is the isotropic pressure, and the PSTF tensor 
TTab = T(^ab) is the anisotropic stress. 

The remaining gauge-invariant variables are formed from the Wcyl tensor Cabcd which vanishes in an exact FRW 
universe because these models are conformally flat. The ten degrees of freedom in the 4-dimensional Weyl tensor can 
be encoded in two PSTF tensors: the electric and magnetic parts defined respectively as 

Eab = Cabcdu'^u'' , (II.8) 
Hab = \Cacstu'v"*bdu'' , (II.9) 

where rjabcd is the 4-dimensional covariant permutation tensor. 



III. FIELD EQUATIONS OF THE BRANEWORLD 

In a recent paper, Maartens introduced a formalism for describing the non-linear, intrinsic dynamics of the 
brane in Randall-Sundrum type II braneworld models in the form of bulk corrections to the 1-1-3 covariant propagation 
and constraint equations of general relativity. This approach is well suited to identifying the geometric and physical 
properties which determine homogeneity and anisotropy on the brane, and serves as a basis for developing a gauge- 
invariant description of cosmological perturbations in these models. 

An important distinction between braneworlds and general relativity is that the set of 1+3 dynamical equations 
does not close on the brane. This is because there is no propagation equation for the non-local effective anistropic 
stress that arises from projecting the bulk Weyl tensor onto the brane. The physical implication is that the initial value 
problem cannot be solved by brane-bound observers. The non-local Weyl variables enter crucially into the dynamics 
(for example, the Raychaudhuri equation) of the intrinsic geometry of the brane. Consequently, the existence of these 
non-local effects leads to the violation of several important results in theoretical cosmology, such as the connection 
between isotropy of the CMB and the Robertson- Walker geometry. 

The field equations induced on the brane are derived by Shiromizu et al [ p6[ using the Gauss-Codazzi equations, 
together with the Israel junction conditions and Z2 symmetry. The standard Einstein equation is modified with new 
terms carrying the bulk effects on the brane: 

Gab = -Agah + K^Tab + K.'^Sab " £ab , (ULl) 
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where = iir/Mp^ . The energy scales are related to each other via 



A 



1 
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(IIL2) 
(III.3) 



where A is the cosmological constant in the bulk and A is the tension of the brane. The bulk corrections to the 
Einstein equations on the brane are made up of two parts: (i) the matter fields which contribute local quadratic 
energy-momentum corrections via the symmetric tensor Sab] and (ii) the non-local effects from the free gravitational 
field in the bulk transmitted by the (symmetric) projection Eab of the bulk Weyl tensor. The matter corrections are 
given by 



Sab — Y^Tc'^Tab — -TacT'^b + ^5ah [ST'crfT'^'' — {Tc'^)^] . 



(IIL4) 



We note that the local part of the bulk gravitational field is the five dimensional Einstein tensor Gab, which is 
determined by the bulk field equations. Consequently, Sab transmits non-local gravitational degrees of freedom from 
the bulk to the brane that includes both tidal (or Coulomb), gravito- magnetic, and transverse traceless (gravitational 
wave) effects. 

The bulk corrections can all be consolidated into an effective total energy density, pressure, anisotropic stress and 
energy flux. The modified Einstein equations take the standard Einstein form with a re-defined energy-momentum 
tensor: 



Gab 



^Agab + , 



(111.5) 



where 



rptot rp 

^ab — J-a, 



T^ab 



jSab 



(111.6) 



Decomposing Sab irreducibly with respect to u"- by analogy with Eq. (ttlj) |21, E2l p3| 



Sab 



K 



UUaUb + 2u^aQb) + ^/^a6 + "Pab 



(III.7) 



(the prefactor is included to make e.g. U have dimensions of energy density), it follows that the total density, pressure, 
energy flux and anisotropic pressure are given as follows: 



P 

ptot 



24 



(2p2-3^'^Vb)+Z^ 



— (4pg„-6^afc(z'') + Qa 



— {-(p + 'iP)-Kab - iT^cl.aT^b)" + M{a(lb) } + V ab 



(III.8) 

(III.9) 
(Ill.fO) 
(Ill.fl) 



For the braneworld case, it is useful to introduce an additional dimensionless gradient which describes inhomogeneity 
in the non-local energy density lA : 



Ta = -DaU 
P 



(III.f2) 



The Gauss-Codazzi scalar equation for the 3-curvature defined by TZ is given by 

4 

K " 

,rp- + 2\ - 
6 V K 



2A 



(III.f3) 
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where 

7^EE (3)^^ = /j^b (3)^^^^ (III. 14) 

with ^^■'-Rafc the intrinsic curvature of the surfaces orthogonal to u"'^. In FRW models the Gauss-Codazzi constraint 
reduces to the modified Friedmann equation 

where the 3-curvature scalar is 7?. = 6K/a^. In non-flat models {K 7^ 0) 7?. is not gauge-invariant since it does not 
vanish in the FRW limit. However, the comoving projected gradient 

Vt ^ ^D,TZ (III.16) 

is a gauge-invariant measure of inhomogeneity in the intrinsic three curvature of the hypcrsurfaces orthogonal to u"". 



IV. LINEARISED SCALAR PERTURBATION EQUATIONS FOR THE TOTAL MATTER VARIABLES 

A. Local and non-local conservation equations 

Based on the form of the bulk energy-momentum tensor and Z2 symmetry, the brane energy-momentum tensor is 
still covariantly conserved: 

v^r,fc = . (IV. 1) 

The contracted Bianchi identities on the brane ensure conservation of the total energy-momentum tensor, which 
combined with conservation of the matter tensor gives 

VSab = K^VSab ■ (IV.2) 

The longitudinal part of Sab is sourced by quadratic energy-momentum terms including spatial gradients and time 
derivatives. As a result any evolution and inhomogeneity in the matter fields would generate non- local Coulomb-like 
gravitational effects in the bulk which back react on the brane. The conservation equation (IV.l) implies evolution 
equations for the energy and momentum densities, and these are unchanged from their general relativistic form. To 
linear order in an almost-FRW brane cosmology we have 

p + Q{p + P)+D''qa = , (IV.3) 

and 

qa + ^eg„ + (p + P)Aa + DaP + D^ab - . (IV.4) 



and 



The linearised propagation equations for U and Q follow from Eq. ( [V.2 ) (see Ref. |22||): 

it + ^eU + D^Qa^O , (IV.5) 

Qa + l^Qa + \DaU + O'^Vab + \u = ^(p + P) {-2DaP + SD^^b + ^Qq^) • (IV.6) 

Taking the projected derivative of Eq. ( [IV. 3 ) we obtain the propagation equation for at linear order: 

pKa + (p + P){Za + aQAa) + aDaD\b + aSDaP - ePAa = 0. (IV.7) 



If the vorticity is non- vanishing flow-orthogonal hypcrsurfaces wiU not exist, and Ti cannot be interpreted as the spatial curvature scalar. 
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From equation ([V.5), we obtain the evolution equation of the spatial gradient of the non-local energy density: 

T„ = ( - - ^ I eT„ - ^-(Z„ + aQAa) - -DaD^Qt. (IV.8) 



P 1 

7^3 



4:U 

3 p 



From the propagation equations for U and Q it can be seen that the energy of the projected Weyl fluid is con served 
while the momentum is not conserved; rather it is driven by the matter source terms on the right of Eq. (IV. 2). Note 
that no propagation equation for Vat is implied so the set of equations will not close. 



B. Propagation and constraint equations 

In this section we give the linearised gravito-magnetic and gravito-electric propagation and constraint scalar equa- 
tions on the brane, which follow from the Bianchi identities, and the equations for the kinematic variables (Jab, and 
Q and its gradient 2a which follow from the Ricci identity. 

For scalar perturbations, the magnetic part of the Weyl tensor Hab and the vorticity tensor ujab vanish identically. 
The electric part of the Weyl tensor Eat and the shear (Tab need not vanish. The non-vanishing variables satisfy the 
following propagation and constraint equations on the brane: 

1. Gravito-electric propagation: 

Eab + QEab + ^l^^iP + P)<^ab + ^K^D^aqb) + '^n'^OT^ab + ^H^T^ab 



1 fk 
72 [k 



[n'^i-Gpip + P)<Jab + 3(p + 3P)7Tab + 3(p + 3P)7T,b 



GpD^aqb) + e[p + 3P]TTab} ~ ^SUffab - SGVab " SGD^aQb) " l^QVab] ; 



2. Shear propagation: 



(IV.9) 



^ab + ^e<Tab + Eab - \>^^^ab - D/^aAb) = ^ ) l^-iP + SP)nab] + l2Vab} \ (IV.IO) 



3. Shear constraint: 



2 1 f 

D^'aab - ^DaQ + K'qa = {k^ pqa + QQa) \ 



(IV. 11) 



4. Gravito-electric divergence: 

1 



D'^Eab + -i^^D'nab 



, 1 2^ 1 f a 

-K-DaP + o = 7^ - 

3 3 48 V K 



-pQqa + 2{p + iP)D\ab + -^pDaP 



IGDaU -16eQa-2AD''Vab 



(IV. 12) 



5. Modified Raychaudhuri equation 

e 



1 r-K^' 



ie2 _^^^^p + sp) + A-^i^] [«V(2p + 3P) + 12U] + D^Aa ; 



(IV. 13) 



6. Propagation equation for the comoving expansion gradient Za which follows from Eq. ( 1V.13| ) 



1 fk-^' 



Za + |eZa - aQAa + '^aDa{p + 3P) - aDaD'Ab = ( ^ ) {^i^aDa[p{2p + 3P)] + l2aDaU} . (IV.14) 
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The spatial gradie nt of th e 3-cu rvatur e scalar is an auxiliary variable. It can be related to the other gauge-invariant 
variables using Eqs. ( III.13| ) and ( III. 16 ): 



6 



- \ 4 
K 



(IV. 15) 



Taking the time derivative of Eq. (IV. 15), commuting the spatial and temporal derivatives, and then making use 
of Eqs. ( IV.13| ) and (IV. 14), we obtain the evolution of the spatial gradient of the 3-curvature scalar: 



1 



Va + ^QVa + ^niZa + aBAa) + -BaDaD'Ab = -{k' + -k^p ) aDaD'qb -[-] aDaD'Q 



1 



~ \ 4 
K 



(IV. 16) 



In general relativity, propagating rja is a useful device to avoid numerical instability problems when solving for 
isocurvature modes in a zero acceleration frame (such as the rest-frame of the CDM) p7| . 



V. COSMOLOGICAL SCALAR PERTURBATIONS IN THE BRANEWORLD 

The tensor-valued, partial differential equations presented in the earlier sections can be reduced to scalar-valued, 
ordinary differential equations by expanding in an appropriate complete set of eigentensors. For scalar perturbations 
all gauge-invariant tensors can be constructed from derivatives of scalar functions. It is thus natural to expand in 
STF tensors derived from the scalar eigenfunctions Q^*^^ of the projected Laplacian: 

satisfying (j^*^^ — 0(1)^- We adopt the following harmonic expansions of the gauge-invariant variables: 





k 


Za 








Eab - 


k 


(Jab 


k 




k 


Aa 


k " 



(V.2) 



Here Va is the 3-velocity of species i relative to for the CDM model considered here we shall make use of Va 
for baryons b and CDM c. For photons 7 and neutrinos v we continue to work with the momentum densities which 
are related to the peculiar velocity of the energy frame for that species by e.g. q^"^ — {4:/3)p'^'^'^v^'' in linear theory. 
The scalar expansion coefficients, such as A^'-* are first-order gauge-invariant variables satisfying e.g. D'^A^j!^ — 0{2). 
Projected vectors q'^'' and STF tensors Q^^^ are defined by p^ : 

2 (V.3) 

so that 



^ The notation 0{n) is short for 0(e") where e is some dimensionless quantity characterising the departure from FRW symmetry. 
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We expand the non-local perturbation variables in scalar harmonics in the following manner: 

k 

Qa^Y.PQ^Qa\ (V.5) 

k 

k 



In addition, we can expand the projected gradient of the 3-curvature term: 



fe 



The form of this expansion is chosen so that if wc adopt the energy frame (where qa — 0) the variable "qk coincides 
with the curvature perturbation usually employed in gauge-invariant calculations. 



A. Scalar equations on the brane 



It is now straightforward to expand the 1+3 covariant propagation and constraint equations in scalar harmonics. 
We shall consider the CDM model where the particle species are baryons (including electrons), which we model as an 
ideal fluid with pressure p^^^ and peculiar velocity vli \ cold dark matter, which has vanishing pressure and peculiar 
velocity v'f'\ and photons and (massless) neutrinos which require a kinetic theory description. We neglect photon 
polarization, although this can easily be included in the 1-1-3 covariant framework |29 . Also, we assume that the 
entropy perturbations are negligible for the baryons, so that DaP^^^ = c^DaP^''^ where cf^is the adiabatic sound speed. 
A complete set of 1-1-3 perturbation equations for the general relativistic model can be found in |3^. We extend these 
equations to brancworld models here. 

In the following, perturbations in the total matter variables are related to those in the individual components by 



k ' 



pqk 



RT^k 



(V.7) 



where q^^^ = (1 + P^^^ / P^^^)^^ ^ (iC^ = and tt^*^ vanishes for baryons and CDM. Similarly, the total density and 
pressure are obtained by summing over components, e.g. P = ^^P^^^- It is also convenient to write P = (7 — l)p, 
but 7 should not be assumed constant (in space or time). 
We begin with the equation for the gravito-electric field: 



1 fk 
72 \k 

'k 



Ik 1 1 

--K^p{-fak - qk) + -K^pQ{l - 37)7rfc -I- -K^p-kk 
2 a 2 



P^ilcTk - qk) - 3(p + iP)pnk - 3(37 - 2)p(p7Tfc -t- pTTfc) - (37 - 2)p^e7rfc 



(V.8) 
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{m<Jk - ipQk) - MpVk + pVk) - upevk 



We have written this equation in such a form that every term is manifestly frame-independent. The shear propagation 
equation is 



k ( ■ lr^ 
a \ 3 



- (*fe + Ak) - —piTk = 7^ - [-(37 - 2)K^p^TTk + \2pVk] . 

a I 2 24 \ K / 



(V.9) 



The shear constraint is given by 



K'^PQk 



2 fk 



3 V a 



3if 



1 ^ 



6 V K 



(kV gfe + 6pQfe). 



(V.IO) 
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The gravito-electric divergence is 



1 fk-^ 



TTfc 



16 Vk 



'3p'e(Zfe + -(37-2)p2(l 



+ 16-pTfe - IGOpQk - 16p 
a 



1 - 



fc2 



K^QpQk 

8k 



3K\ k 

a 3 a 



(V.ll) 



The propagation equation for the comoving expansion gradient Za is given by 

2(p(-) Ai^^) + pC^^Ai^^) + (1 + Sc^jpC) Ai'') + p(^) Ai^) 
■ 1^ i {k^[(2p + 3P)pA, + p(3p(^)Ai-^^ + 3pMAM + (2 + 3cf)pW A^ + 2pW A^^')] + 12pT,} . 

The non-local evolution equations for and Qk are 



2k + gOZfc - -eAfc + -Ak + —-j^ 

1 / - \ 4 

1 I k\ a 



(V.12) 



Tfc = ^(37-4)eT,-^e^Afe 



and 

&-i(37-4)eQ. + i^ 



Tfc + 2 ( 1 



3ii: 



AkU 
3 a p' 



AUk„ k^ 

o -2fc + -Sfe, 

6 p a a 



IP { Qqk 



3K 



TTfc — Afe 



The spatial gradient of the 3-curvature scalar is 

2 



1 



3K 



9 -49 
K^P . u p . 

Vk = -^'^k + -J^^k 



1 /k 



Ik 



2UJ '^'-3a^^'' 



(V.13) 



(V.14) 



(V.15) 



and it evolves according to 



3K 



r]k - gQ^fe 



K 



1 /k 



^,2k - 2« Pik ^2 V K 



{K^p^qk + 6pQfe). 



(V.16) 



The evolution equations for the scalar harmonic components of the comoving, fractional density gradients for 
photons, neutrinos, baryons and cold dark matter (CDM) are 



Ai^^ 



-0Afe (photons) , 



k fA 



-QAk (neutrinos) 



P 

-^{Zk-vi'^)-eAk (CDM). 



1+— ; 



The evolution equations for the momentum densities and peculiar velocities are 



•(7) 

9fe = 


3 a 


Qk = 


Ik 
3 


(pW+pW)4'') = 




.(c) 


--Gv 
3 



A 



(7) 

k 



■ 4Afc + 2 1 



3ii:^ 









(b) (7) 



(1-3^64") + 



- -(1 + c^)A^''^ - n ( -v^''^ - 0^^^ 



(V.17) 
(V.18) 
(V.19) 
(V.20) 

(V.21) 
(V.22) 
(V.23) 
(V.24) 
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where the Thomson scattering terms involving the electron density rig and Thomson cross section try arise from the 
interaction between photons and the tightly-coupled baryon/electron fluid. The remaining equations are the propa- 
gation equations for the anisotropic stresses of photons and neutrinos, and the higher moments of their distribution 
functions. These equations can be found in ||30| , and with polarization included in |29| , since they are unchanged from 
general relativity. However, we shall not require these additional equations at the level of approximation we make in 
our subsequent calculations. 



VI. PERTURBATION DYNAMICS IN THE CDM FRAME 



In this section we specialize our equations to FRW backgrounds that are spatially flat'^ and we ignore the effects 
of the cosmological constant in the early radiation-dominated universe. To solve the equations it is essential to make 
a choice of frame u". In Ref. two of the present authors adopted a frame comoving with the CDM. Since the 
CDM is pressure free, this u"" is geodesic {Aa — 0) which simplifies the equations considerably. We shall adopt this 
frame choice here also, though we note it may be preferable to use a frame more closely tied to the dominant matter 
component over the epoch of interest. This can be easily accomplished by adopting the energy frame {qa = 0). For 
completeness, we give equations in the energy frame in the appendix. 

We neglect baryon pressure (c^ and P*^**^ 0) and work to lowest order in the tight-coupling approximation 
(neCTT oo; see e.g. Ref. |pl[). At this order the energy frame of the photons coincides with the rest frame of the 
baryons, so that v^^ ~ 3gi^(4p^^)), and all moments of the photon distribution are vanishingly small beyond the 
dipole. 

With these approximations and frame choice we obtain the following equations for the density perturbations of 
each component: 

Ai^^ = -^Qz.-gI^^) (photons), (VI.l) 
Ai:^' = (^^^ - 1k^) (neutrinos) , (VI.2) 
Ai') =-^(Zfe-t;f)) (baryons), (VI.3) 
A(^^)=-^Zfe (CDM). (VI.4) 
The equations for the peculiar velocities and momentum densities are 

(4p(^) + 3p(''))g[^) = -l^p(^)A(^) - p^'^eql?'^ , (VI.5) 
^-fAr+2.r) , (VI.6) 



3 a 

„(=) _ n Jb) _ „ (7) 



along with w^, = and vj, — 3g^ /4. The latter equation, together with Eqs. (VI.l) and (VI.3), implies that 



^if^ = 3A^''''/4 so that any entropy perturbation between the photons and baryons is conserved while tight coupling 



holds. The effects of baryon inertia appear in Eq. (VI.5) because of the tight coupling between the baryons and 
photons. 

The constraint equations are found to be: 

{Zk-<Jk)-U^) («W + 6pQfe) = 0, (VI.7) 



^ More generally, curvature effects can be ignored for modes with wavelength much shorter than the curvature scale, k S> y^j-R"!, provided 
the curvature does not dominate the background dynamics. 
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and 



$fe - K> ( - ) (Afc - TTfc) + K^Qpqk 



1 /k-^' 



16 V K 



-^Qphk + ^^p' [(37 - 2)7rfe + 2Ak] 



(VI.8) 



16 - p(Tfc - T'fc) - IGOpQ, 



The propagation equation for the comoving expansion gradient in the CDM frame is 



1 a 

1 / K \ a 

^ "l2 



{n\2p + iP)pAu + p(3p(-)Ai-) + Sp^Ai'') + (2 + 3c^)p(^)Ai^) + 2p(^)Ai^))] + 12pT,} 



(VI.9) 



The variables $fc and (7^ can be determined from the constraint equations so their p ropa gation equations are not 
independent of the above set. The propagation equation for (f>fc is unchanged from Eq. ( V.8 ) since that equation was 
already written in frame-invariant form. The propagation equation for the shear in the CDM frame is 



1 fk^' 



h(37-2)KV 7rfc + 12p7'fc] . 



(VI. 10) 



Finally we have the non-local evolution equations for and Qk which in the CDM frame become 



tfc-i(37-4)eT,.-^--Zfe + -Qfc 
6 Spaa 



(vi.ii) 



and 



Qk - i(37 - 4)eQfe + \-{T, + 2Vk) = ^7P 
3 6 a b 



eqfc + -(TTfc- Afc ; 

a 



(VI. 12) 



VII. SOLUTIONS IN THE RADIATION-DOMINATED ERA 



We now use the above equations to extract the mode solutions of the scalar perturbation equations in the radiation- 
dominated era, 7 = 4/3. To simplify matters, as well as neglecting the contribution of the baryons and CDM to the 
background dynamics, we shall only consider those modes for which Dap'^''^ and DaP^'^^ make a negligible contribution 
to the total matter perturbation Dap. This approximation allows us to write the total matter perturbations in the 
form 

(p(^) + p('^))Afc = p(^)Ai^) + p('^^Ai''\ (p(^) + /'^))gfc = p^^ki'^ + P^^ki^K (VII.l) 

and effectively removes the back-reaction of the baryon and CDM perturbations on the perturbations of the spacetime 
geometry. We note that in making this approximation we lose two modes corresponding to the baryon and CDM 
isocurvature (density) modes of general relativity, in which the sub-dominant matter components make significant 
contributions to the total fractional density perturbation (which vanishes as t — *■ 0). However, for our purposes the 
loss of generality is not that important, while the simplifications resulting from decoupling the baryon and photon 
perturbations are considerable. We also neglect moments of the neutrino distribution function above the dipole (so 
there is no matter anisotropic stress). This approximation is good for super- Hubble modes, but fails due to neutrino 
free streaming on sub-Hubble scales. 

We shall also assume that the non-local energy density U vanishes in the background for all energy regimes pl| . 
Physically, vanishing U corresponds to the background bulk being conformally flat and strictly Anti-de Sitter. Note 
that U — in the background need not imply that the fluctuations in the non-local energy density are zero, i.e. 

Ta^O. 
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With the above conditions the following set of equations are obtained: 



K^p / k 



k2 a 



-]Qk + iHVk - Vk 
a 



i&k + HcTk) 



6 



k2 a 









3 a 


= 0, 








3a 


= 0, 




k 
h - 

a 






= 0, 




k 










h - 
a 






= 0, 



where recall H = 0/3. For the constraint equations we find 



A 



18 p 



k2 A' 



3Hqk 



6 p 
^A 



(Tfc - Vk) - 3HQk 



Finally the non-local evolution equations are found to be : 



Tfc = -Qk , 
a 



9Qfe + 3 - (Tfe + 2Vk) = -2kV -Afe - 3Hqk 



It is easy to show by propagating the constraint equations that the above set of equations are consistent. 
By inspection, there is a solution of these equations with 



= 0, 



2k ~ 



3H 
A 

Cfc = , 

o 

3 a 
AA 

3 a ' 



- 1 



A 

a 



it) 
Ik 

it' 
Tfc 



A 



0, 



Qfc = 0, 
Vk^Q, 



(VII.2) 

(VIL3) 

(VII.4) 

(VIL5) 
(VIL6) 
(VIL7) 

(VIL8) 

(VII.9) 
(VII. 10) 

(VII. 11) 
(VII.12) 

(VII. 13) 
(VII. 14) 

(VII. 15) 
(VII. 16) 
(VII. 17) 
(VII. 18) 

(VII. 19) 

(VII.20) 
(VII.21) 
(VII.22) 



where A is a constant. This solution describes a radiation-dominated universe that is exactly FRW except that the 
CDM has a peculiar velocity v'^a'^ = {A/a)Q^a'^ relative to the velocity of the FRW fundamental observers. [This form 
for v'^a'^ clearly satisfies Eq. (V.24) with Aa = 0.] Such a solution is possible since we have neglected the gravitational 
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effect of the CDM (and baryon) perturbations in making the approximations in Eq. (VII.l). The same solution arises 
in general relativity Including the back-reaction of the CDM perturbations, we would find additional small 

peculiar velocities in the dominant matter components which compensate the CDM flux. We shall not consider this 
irregular CDM isocurvature velocity mode any further here. 

Another pair of solutions are easily found by decoupling the photon/neutrino entropy perturbations. Introducing 
the photon/neutrino entropy perturbation (up to constant) A2 and relative flux q2: 



the equations for A2 and q2 decouple to give 



A2 

92 



(f) 



(VII.23) 



A, 



-92 



92 



Ik 

L 

3 a 



0, 
. 



Switching to conformal time {dr — dt/a) we can solve for A2 and q2 to find 

' kr^ 



92 (r) 
A2(t) 



B cos 



B sin 



fcr 

T 

fcr 

Y 



- Csin 
Ccos 



(VII.24) 
(VII.25) 

(VII.26) 
(VII.27) 



The constants B and C label the neutrino velocity and density isocurvature modes respectively [^0[ in which 
the neutrinos and photons initially have mutually compensating peculiar velocities and density perturbations. The 
perfect decoupling of these isocurvature modes is a consequence of our neglecting anisotropic stresses (and higher 
moments of the distribution functions) and baryon inertia. 

Having decoupled the entropy perturbations, we write the remaining equations in terms of the total variables A^ 
and qk- The propagation equations for the non-local variables and Qk are redundant since these variables are 
determined by the constraint equations ( VII.9D and (VII.IC): 



k2 A 



3 _ 1 



$fc + 2H {^^ [Zk - <Jk) - (1 + Afc + ^^Vk. 



(Zfe — cTfc) — 



2 



9fc- 



(VII. 28) 
(VII.29) 



Substituting these expressions in the right-hand sides of Eqs. (VII. 2) and (VII. 2) we find 



H 



A, 



1 

" 3 
2H 



{Zk - CTfc) 

(Zfe - o-fc) 

- 1 ^fe 



k /4 



Ak + -[ -Z. 

qk 4 



9fc 



1 k 

3 a 



A, = 



Ki A 


(VII.30) 


K^A 


(VII.31) 


3 p 

^A ' ' 


(VII.32) 


0, 


(VII.33) 


. 


(VII.34) 



These equations describe the evolution of the intrinsic perturb ations to the brane. The usual general relativistic 
constraint equations are now replaced by the constraints (VII. 28) and ( VII.29| ) which determine two of the non-local 
variables. The lack of a propagation equation for Vk reflects the incompleteness of the 1-1-3 dimensional description 
of braneworld dynamics. 

In the following it will prove convenient to adopt the dimensionless independent variable 



k 

Ha ' 



(VII.35) 



14 



which is (to within a factor of 27r) the ratio of the Hubble length to the wavelength of the perturbations. Using the 
(modified) Friedmann equations for the background in radiation domination, and with U = 0, we find that 



dx k f 2 



3p/A \ 



dt a\2 + p/\ J 



(VII. 36) 



The relative importance of the local (quadratic) braneworld corrections to the Einstein equation depends on the 
dimensionless ratio p/A. In the low-energy limit, p <C A, the quadratic local corrections can be neglected although 
the non-local corrections Eab may still be important. In the opposite (high-energy) limit the quadratic corrections 
dominate over the terms that are linear in the energy- momentum tensor. We now consider these two limits separately. 



A. Low-energy regime 



In the low-energy regime we have dx/dt fa k/a and x « fcr. The total energy density p is proportional to x~'^. 
Denoting derivat ives wit h respect to x with a prime, using p <^ X, and assuming that we can neglect the term involving 
{p/X)Ak in Eq. ( VII. 31 ) compared to the other terms, we find 



3a;2$fc -I- 3x$fc + (6 + x^)ak - x^Z^ 
x^Z'f. + 3xZk + 2a;^$fc - 2x(7k 



27 



x^o-fe + xak + x^^k 
4 



Ai 



3 



1 







V, 



3 



:Ak = 



Combining these equations we find an inhomogeneous, second-order equation for 

3x'S>'l + 12$; + a;$fc = Fkix), 

where 



(VII.37) 
(VII. 38) 
(VII. 39) 
(VII.40) 
(VII.41) 

(VII.42) 

(VII.43) 



In general relativity the same second-order equation holds for but with Fk{x) = 0. 

The presence of terms involving the non-local anisotropic stress on the right-hand side of Eq. (VII.42) ensure that ^k 
cannot be evolved on the brane alone. The resolution of this problem will require c areful a nalysis of the bulk dynamics 
in five dimensions. In this paper our aims are less ambitious; we shall solve Eq. ( VII.42| ) with Vk — 0. Although we 
certainly do not expect Vab = 0"*, the solutions of the homogeneous equation may still prove a use ful star ting point for 
a more complete analysis. For example, they allow one to construct Green's functions for Eq. ( VII.42 ) which could 
be used to assess the impact of specific ansatze for Vab |33l . 



* We have not investigated the consistency of the condition Vab = with the five-dimensional bulk dynamics in the presence of a perturbed 
brane. 
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With Vk = we can solve Eqs. (|VII.37|)-(|VII.4lD analytically to find 



3 



3 sin 



X 



C3[b + X 



Ok 
Zk 

Afc = C4 COS 



—= I — X^/icOS ( —pz I 



C2 COS ( —= I + ci sin ( ) 



£3 

a; 



3 COS f I + a:\/3 £ 



2) , 6V3 



c,cos,^]-c.sin(-^] 



C2 COS 



V3j 
4a/3 2 ^ 

7!. 



• ^ ^ ^ ^ 4 



a; 



C2Cosl^j+c,sin^^ 



a; 



qk 



4c3 4x 

\/3 """^ V\/3y' \/3""V\/3y' ' 3^ 75 



ci cos I — = I — C2 sin 



%/3) 



C2 cos I + ci sm I ^= I 



V3j 



(VII.44) 
(VII.45) 
(VII.46) 

(VII.47) 

(VII.48) 



The mode labelled by C3 is the CDM velocity isocurvature mode discussed earlier. The modes labelled by ci and C2 are 
the same as in general relativity; they describe the adiabatic growing and decaying solutions respectively. However, 
in the low-energy limit we also find two additional isocurvature modes (04 and C5) that are not present in general 
relativity. These arise from the two additional degrees of freedom and Qk present in the braneworld model (with 
Pk = 0). The mode C4 initially has non-zero but compensating gradients in the total matter and non-local densities, 
and C5 initially has compensated energy fluxes. Formally these isocurvature solutions violate the assumption that 
the term involving {p/X)Ak be negligible compared to the other terms in Eq. (VII. 31) since all other terms vanish. 



In practice, there will be some gravitational back-reaction onto the other gauge-invariant variables controlled by the 
dimensionless quantity p/A, but the general character of these isocurvature modes will be preserved for p/ X <C 1. 



B. High-energy regime 

We now turn to the high-energy regime, where the quadratic terms in the stress-energy tensor dominate the (local) 
linear terms. In this limit the scale factor a oc t^/^. The modification to the expansion rate leads to an increase 
in the amplitude of scalar and tensor fluctuations produced during high-energy inflation ||3^ . With U = in the 
background, and p 3> A, the Hubble parameter is approximately 



36 



-4 2 
K p , 



and dx/dt « 3k /a. In terms of conformal time r, a; ~ 3fcr. 



1. Power series solutions for the high-energy regime 



(vn.49) 



It is convenient to rescale the non-local variables by the dimensionless quantity n'^p. Thus we define 



Tk = 



Qk 



Vk 



K'^p 

Qk_ 

K^P 

Vk_ 

K.'^P 



(vn.50) 
(vn.51) 

(VII. 52) 
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The fractional total (effective) density perturbation and energy flux can be written in terms of the barred variables 
[e.g. Tq = Tq/(k^p)] in the high-energy limit as 



aDap" 



2(A, + 6Ta), 



P 



Making the high-energy approximation p ^ A in Eqs. ( VII.30)-( VII.34), we obtain 



9x2$fc + ix^k + (12 + x^)ak - x^Zk 



54 



X " 



2xak + 2x'^^k 



a1- 



12Afc 



1 4^ 
3* + 9^" 



Ik 



1 ^ 
9 



Vk 
18— , 

0, 
. 



The non-local quantities and are determined by the constraints 

18 18 6 



-X (Zfc - Ok] 



(Ik- 



(VII.53) 
(VII.54) 

(VII.55) 
(VII.56) 
(VII.57) 
(VII. 58) 
(VII.59) 

(VII.60) 
(VII.61) 



54 ' " "'6 

We can manipulate Eqs. ( VII. 55| )-( VII.59) to obtain a fourth-order equation for the gravitational potential <I>fc: 



where 



729x- 



54 



5x4 



3888x- 



(1782 + 54x^ 



dx^ 



lUx^ + i90 + x^)^k=Fkix) , 



(VII. 62) 



Ffc(x) = J 243x 



1 9^ 



SlOx-' 



dx^ 



18x^(135 -f2x^ 



9x^ 



30x(162-f x^ 



9x 



x** + 30(162 + x^)]7'fc 

(VIL63) 



Since we do not have an evolution equation for Vk we adopt the strategy taken in the low-energy limit and look 
for solutions of the homogeneous equations {Vk = 0). In principle one can use these solutions to construct formal 
solutions of the inhomogeneous equations with Green's method. 

To solve Eq. (VII. 62) with = we construct a power series solution for <I>fe(x): 



$fc(x) =x™^a„a 



ra=0 



where oq 7^ 0. The indicial equation for m is 

771(171 — l)(3m -|- 5)(3m — 4) = . 



(VII.64) 



(VII.65) 



For each value of m we substitute into Eq. (VII. 62) and solve the resulting recursion relations f or the | a „|. We then 
obtain the other gauge-invariant variables by direct integration. The original set of equations ( VII. 55)-( VII.59) has 
five degrees of freedom, so we expect one additional solution with ^k = 0. This solution is the CDM isocurvature 
solution discussed earlier, and has a finite series expansion: 



= 0, 

Zk = Cx-3(12-fx2) 
Afc = ACx-i , 



qk 



-Cx-^ , 



(VII.66) 
(VII.67) 
(VII.68) 
(VII.69) 

(VII.70) 
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where C is a constant. The non-local variables vanish. 
The first two terms of the mode with m — are 



CTfe 

Qk 



bi 
bi 
bi 
bi 
bo 
bi 



198 

1 1 c 

—X H 

4 396 



1 2 1 4 
-X X 

6 864 
1 . 



1 



162 



972 



38880 
1 

249480" 
1 



bi { x^ + 

' 243 17820 



(VII.71) 
(VII.72) 
(VII.73) 
(VII.74) 
(VII.75) 
(VII.76) 
(VII.77) 



where bi is a constant. The form of this solution is similar to the adiabatic growing mode of general relativity. 
The mode corresponding to m = 1 is 



Qk 
Tfc 
Qk 



b2 
b2 
b2 
b2 
b2 
b2 
b2 



13 
1890 



1890 



-r 

72 12 



7 35 
18 1 , 

--"+15^^ 



6 + -x^ 



"+30"^ 



(VII.78) 
(VII.79) 
(VII.80) 
(VII.81) 
(VII.82) 
(VII.83) 
(VII. 84) 



with 62 a constant. As t ^ there are non-zero but compensating contributions to the effective peculiar velocity 
from the matter and the non-local energy fluxes. The contributions of these components to the fractional 
total density perturbation a-DaP*°'/p*°* vanish as i 0. It follows that this solution describes an isocurvature velocity 
mode where the early time matter and non-local (Weyl) components have equal but opposite peculiar velocities in 
the CDM frame. The existence of such isocurvature modes was predicted in Refs |l^ and for large-scale density 
perturbations. 
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The mode corresponding to to = — | is singular as < ^ (it is a decaying mode): 

5 



<Jk = bsx-i (^l + ^x^^ , (VII.86) 

Z, = 63 f llxt - ^^x^^ , (VII.87) 
^ 1590435 y ' ^ ' 

^^-^^(-4^^ + 43^55^^)' (^"-^^^ 
^^ = ^3(^x^-^xf), (VII.89) 

T. = 63(-^.i + ^.^), (VII.90) 

Qfe = &3 f-— a;3 + — ^ . (VII.91) 
^ V 54 4860 J ^ ' 

A similar mode is found in general relativity but there the decay of <i>fe is more rapid ($fc cx x~^) on large scales. 
Finally, for m = | we have 

$fe = 642;^ [i - -^^^^ ' (vn.92) 

= (VII.93) 

^M^(f-ig.^), (VII.94) 

Afc = 64X^(^-^ + 1x2) , (VII.95) 

Qk = b,x-^ (^y^x - -^x^^ , (VII.96) 

T(x) = bixi Q + ^x^) , (VII.97) 

In this mode the universe asymptotes to an FRW (brane) model in the past as t — > 0. Note that this requires careful 
cancellation between aDaP^°^ / p^°^ and to avoid a singularity in the gravitational potential $fc (which would 
diverge as x~'^^^ without such cancellation). Like the velocity isocurvature mode (to = 1) discussed above, this mode 
has no analogue in general relativity. 

VIII. A COVARIANT EXPRESSION FOR THE TEMPERATURE ANISOTROPY 

In this section we discuss the line of sight solution to the Boltzmann equation for the scalar contribution to the 
gauge- invariant temperature anisotropy Sxie) of the CMB in braneworld models. We employ the 1+3 covariant 



approach, and show that our result is equivalent to that given recently by Langlois et al |11 using the Bardeen 
formalism. 

Over the epoch of interest the individual matter constituents of the universe interact with each other under gravity 
only, except for the photons and baryons (including the electrons), which are also coupled through Thomson scattering. 
The variation of the gauge- invariant temperature perturbation S^ie), where e° is the (projected) photon propagation 
direction, along the line of sight is given by the (linearized) covariant Boltzmann equation (valid for scalar, vector. 
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and tensor modes) 



(VIII.l) 



where the prime denotes the derivative with respect to a parameter A defin ed alon g the hne of sight by dX — —Uadx'^. 

Following the steps in Ref. [Q, we expand the right-hand side of Eq. (VIII.l) in scalar harmonics and integrate 
along the line of sight from the early universe to the observation point R. Neglecting effects due to the finite thickness 
of the last scattering surface, on integrating by parts we find that the temperature anisotropy involves the quantity 



Ik. 



HA, 



f- 



I 



(VIIL2) 



integrated along the line of sight (after multiplying with Q^*^-*). In simplifying Eq. ( VIII.2| ) we have made use of the 



derivative of the shear propa gation equat ion (V.9), substituted for qk and Zk from equations (V.8) and (V.IO), and 



finally used equations ( III. 15 ) and ( IV. 13 ). The quantity / is the total sum of all the braneworld corrections: 



(VIII.3) 



where 



h 
h 



1 fh-' 



24 \k 
1 fk 



72 Vk, 
(37 - 2)p2e7rfc 



[-(37-2)KV'7rfc + 12p7'fe] , 



6 ^-^ ip^dk ~ 3(p + 3P)pTrk - 3(37 - 2)p{p'hk + pT^k) ~^{^] P^l^ 
48 Uuk - 36{pVk + pVk) + 36 pQu - UpOVk 



(VIII.4) 



1 /k 



12 V K 
1 



[(37- 1)kV^ + 12U 



18 o \K J 24 V K ,' 



A lengthy calculation making use of the propagation and constraint equations shows that / = 0. The final result for 
the temperature anisotropics is then 



[STie)]H = -J2 



+ Ak ] Q("-) 



+ 



3 

16 



k k ''^A 



(VIII.5) 



where the event A is the intersection of the null geodesic with the last scattering surface 

In retrospect, one could re-derive the result for the temperature anisotropy in braneworld models much more simply 
by retaining the effective stress-energy variables p'"' 



ptot ^ ^to t g^j^j ^tot jj-^ ^Yi^ propagation and constraint equations 



used in the manipulation of the left-hand side of Eq. ( VIIL2), rather than isolating the braneworld contributions. 

If we adopt the longitudinal gauge, defined by aab = 0, we find that the electric part of the Weyl tensor and the 
acceleration are related by = —Ak if the total anisotropic stress Trjj*^' vanishes. It follows that in this zero shear 
frame we recover the result found by Langlois et al jl^. 

Regarding the imprint of braneworld effects on the CMB, we note several possible sources. Once the universe 
enters the low-energy regime the dy namic s of the perturbations are essentially general relativistic in the absence of 
non-local anisotropic stress (see Sec. VILA). If Tab really were zero, the only imprints of the braneworld on the CMB 
could arise from modifications to the power spectrum (and cross correlations) between the various low-energy modes. 
Since there are two additional isocurvature modes in the low-energy universe due to braneworld effects, it need not 
be the case that adiabatic fluctuations produced during high-energy (single-field) inflation give rise to a low-energy 
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universe dominated by the growing, adiabatic, general-relativistic mode. The possibility of exciting the low-energy 
isocurvature (brane) modes from plausible fluctuations in the high-energy regime is worthy of further investigation. 
In practice we do not expect Vab = 0. In this case the non-local anisotropic stress provides additional driving terms 
to the dynamics of the fluctuations, and we can expect a significant manifestation of five-dimensional Kaluza-Klein 
effects on the CMB anisotropics. 

IX. CONCLUSION 

In this paper we have discussed the dynamics of cosmological scalar perturbations in the braneworld scenario from 
the viewpoint of brane-bound observers making use of the 1-1-3 covariant approach. We only considered matter 
components present in the ACDM model, but it is straightforward to include other components such as hot dark 
matter. 

We presented approximate, analytic solutions for the fluctuations in the low-energy universe under the assumption 
that the non-local anisotropic stress was negligible. We obtained two additional isocurvature modes not present in 
general relativity in which the additional density gradients or peculiar velocities of the total matter are compensated 
by fluctuations in the non-local variables. In practice we do not expect the non-local anisotropic stress Vab necessarily 
to be negligible; in this case our solutions to the homogeneous equations should form a useful starting point for 
the construction of solutions to the driven equations. By adopting a four-dimensional approach our presentation is 
necessarily limited. In particular, we cannot predict the evolution of Vat on the brane. However, the four-dimensional 
approach should be well-suited to a phcnomcnological description of these flve dimensional Kaluza-Klcin modes. A 
simple possibility is to adopt an ansatz for the evolution of Vab ||3^ , and this will be explored further in a future 
paper. 

We also presented solutions to the perturbation equations in the high-energy regime where braneworld effects 
dominate. In this limit the gravitational potential satisfles a fourth-order differential equation which we were unable 
to solve analytically (even with Vab — 0). Instead we constructed power series solutions for the case where the non- 
local anisotropic stresses vanish; these should prove useful for setting initial conditions in the high-energy regime when 
performing a numerical solution of the perturbation equations. We found two additional modes over those present in 
general relativity, one of which can be described as a (brane) isocurvature velocity mode. We also showed that the 
adiabatic decaying mode varies less rapidly than in general relativity on large scales. 

The detailed calculation of braneworld imprints on the CMB (in the phenomenological approach discussed above) 
will be described in a future paper. Here we showed with the 1-1-3 covariant approach that the line of sight integral 
for the CMB temperature anisotropics is unchanged in form from general relativity. We also noted that excitation 
of the additional isocurvature modes present in the low-energy universe could provide an additional imprint in the 
CMB, over and above that due to the non-local anisotropic stress. 
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APPENDIX A: ENERGY FRAME EQUATIONS IN THE RADIATION-DOMINATED ERA 

In this appendix we present a complete set of evolution equations for the total matter variables in the matter 
energy frame, Qa = 0. Note that the four-velocity of the energy frame is not necessarily a timelike eigenvector of 
the Einstein tensor in the presence of the non-local braneworld corrections to the effective stress-energy tensor. We 
assume that the matter is radiation dominated, the non-local energy density vanishes in the background, and we 
ignore local anisotropic stresses. We also assume that the baryons and CDM make a negligible contribution to the 
fractional gradient in the total matter energy density and to the energy flux, thus excluding the CDM and baryon 
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isocurvature modes. We also give the evolution equations for the non-local density gradient and energy flux in the 
matter energy frame. 

Denoting the variables in the energy frame by an overbar^, the relevant equations for scalar perturbations are 

L = ^QAa-^Za, (A.Al) 

Ta = --D^Qa , (A.A3) 

p 

Qa = ~t^Qa~^fa~^^Aa~D''Vab. (A.A4) 

6 6a ya 



Solutions of these equations are related to those in the CDM frame (Sec. VII) by linearising the frame transforma- 
tions given in Ref. |^. If the CDM projected velocity is ui^"* in the energy frame, the variables in the CDM frame 
are given by 

' ' (A.A5) 
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r . (A.A6) 

(A.A7) 
(A.A8) 

(A.A9) 

where we have used lA = Q in the background. The CDM peculiar velocity evolves in the energy frame according to 

v'^:^ = + ^Aa. (A.AIO) 
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